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Abstract 



A search for the light pseudoscalar sgoldstino production in the three-body K~ decay K~ — > 7r~7r° P 
has been performed with the "ISTRA+" detector exposed to the 25 GeV/c negative secondary beam of 
the U-70 proton synchrotron. No signal is observed. An upper limit for the branching ratio Br{K' — > 
7r~7r°P), at 90% confidence level, is determined to be ~ 9 ■ 10"^ in the effective mass nip range 
from to 200 MeV/c^, excluding the region near uIt^o where it degrades to ~ 3.5 • 10~^. 
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1 Introduction 



In models with spontaneous supersymmetry breaking the 
superpartners of a Goldstone fermion, pseudoscalar P and 
scalar S goldstinos, should exist. In some versions of 
gravity-mediated and gauge-mediated theories ( for a re- 
cent review see [T] ) one or both of these weakly interact- 
ing bosons (sgoldstinos) are light enough to be observed 
in kaon decays. Moreover, if sgoldstino interactions with 
quarks conserve parity, as in left-right extensions of MSSM, 
and P is lighter than S, so that ms > {niK — rn^) and 
nip < (rriK — 2m^), sgoldstinos can be observed in the de- 
cay K 7r7rP(see Fig.l), rather than in the much better 
constrained K — > irS. The phenomenology of light sgold- 
stinos in this scenario is considered in detail in ||^|. Under 

the assumption that sgoldstino interactions with quarks and gluons violate quark flavor and conserve 
parity, low energy interactions of pseudoscalar sgoldstino P with quarks are described by the La- 
grangian: 

(1) 




Figure 1: 
and pions. 



K decay into sgoldstino 



L = -P ■ (/ig ■ dii-f^dj + h'r'. ■ Uii-f^Uj 



where 



di = {d, s, h) 



(u,c,t) 



and with coupling constants hij proportional to the left-right soft terms in the mass matrix of squarks: 
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where the scale of supersymmetry breaking is denoted as ^/F. 

The 90% confidence level (CL) constraints on the flavor-violating coupling of sgoldstinos to quarks 
evaluated using the K°^ — Kg mass difference and CP violating parameter e in the neutral kaon system 
are: j/ifal < 7 ■ 10"^; \KQ{h^2)^m{h^^)\ < 1.5 • 10"^'^. It has been shown O that, depending on 
the phase of sgoldstino-quark coupling, these constraints result in the following 90% CL upper limits 
on the branching ratio: Br(K^ n'lr^P) < 1.5 ■ 10^^ — 4 • lO^'', where the less strong 
limit corresponds to the case of pure real or pure imaginary /if^ . A search for P in charged kaon 
decays is of particular interest when Re(/if^) ~ 0, and when related branching ratio for the decay 
Kl n^n^P is small. 

Light sgoldstinos decay into two photons or into a pair of charged leptons. The two photon decay 
dominates almost everywhere in the parameter space. Depending on the parameter = ^^^^^ 
where M^^ is the photino mass, sgoldstinos will have a very different lifetime. In the present search 
we assume that the sgoldstino is sufficiently long lived to decay outside the detector, i.e. the sgoldstino 
is "invisible". The existing 90% CL upper limit on the branching ratio Br{K^ — > n^-n^P) is 4 ■ 10~^ 
Q, whereas the 90% CL upper limits for the scalar sgoldstino S estimated from the studies of the 
^Ti+vV vary between 0.4 ■ 10^^° and 1.0 • 10^^° in the mass interval from to 110 MeV/c^, 
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Figure 2: The elevation view of the ISTRA+ setup. Mi and M2 are magnets, Cj — Cerenkov counters, 
BPCj — beam proportional chambers, PCj — spectrometer proportional chambers, SP,; — lead glass 
calorimeters, DCj — drift chambers, DTj — drift tubes, HC — a hadron calorimeter. Si — trigger 
scintillation counters, MH — a scintillating hodoscope and MuH — a scintillating muon hodoscope. 



see Fig. 3 in [l6l). The recent upper limit for the decay n^X°, where X° is a neutral weakly 

interacting massless particle, is 0.73 • 10~^°(90% CL)|I3. 

The aim of our present study is to search for invisible pseudoscalar sgoldstinos in the K~ decays. 
The experimental setup and event selection are described in section 2, the results of the analysis are 
presented in section 3, the systematic uncertainties are discussed in section 4 and the conclusions are 
given in the last section. 



2 Experimental setup and event selection 

The experiment is performed at the IHEP 70 GeV proton synchrotron U-70. The ISTRA+ spectrom- 
eter has been described in some detail in our recent papers on K^s fl], K^^, |l8| and n^n°n° decays 
||9ll . Here we recall briefly the characteristics relevant to our analysis. The ISTRA+ setup is located 
in a negative unseparated secondary beam line 4A of U-70. The beam momentum is ~ 25 GeV/c 
with Ap/p ~ 1.5%. The admixture of in the beam is ~ 3%, the beam intensity is ~ 3 ■ 10^ 
per 1.9 sec U-70 spill. A schematic view of the ISTRA-i- setup is shown in Fig. 2. The beam particles 
are deflected by the magnet Mi and are measured by four proportional chambers BPCi — BPC4 with 
1 mm wire spacing, the kaon identification is done by three threshold Cerenkov counters Co — C2. 
The 9 meter long vacuum decay volume is surrounded by eight lead glass rings used to veto low 
energy photons. The 72-cell lead-glass calorimeter SP2 plays the same role. The decay products are 
deflected in the magnet M2 with 1 Tm field integral and are measured with 2 mm step proportional 
chambers PCi — PC3, with 1 cm cell drift chambers DCi — DC3 and, finally, with 2 cm diameter drift 
tubes DTi — DT4. The wide aperture threshold Cerenkov counters C3 ,64 , filled with He, serve to 
trigger electrons and are not used in the present measurement. SPi is a 576-cell lead-glass calorime- 
ter, followed by HC, a scintillator- iron sampling hadron calorimeter. MH is a 1 1x1 1 cell scintillating 
hodoscope, used to improve the time resolution of the tracking system, MuH is a 7x7 cell muon 
hodoscope. 

The trigger is provided by scintillation counters Si — S5, beam Cerenkov counters and by the 
analog sum of amplitudes from last dynodes of the SPi : T=Si ■ S2 ■ S3 ■ S4 ■ Ci ■ C2 ■ C3 ■ S5 ■ S(SPi), 
here S4 is a scintillation counter with a hole to suppress the beam halo, S5 is a counter downstream of 
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the setup at the beam focus, S(SPi) requires that the analog sum to be larger than the MIP signal. 

During the first run in March- April 2001, 363 million of trigger events were recorded. During 
the second physics run in November-December 2001 350 million trigger events were collected with 
higher beam intensity and stronger trigger requirements. This information is complemented by about 
300 M Monte Carlo (MC) events generated using Geant3 [10| for the dominant decay modes. 
Signal efficiency for possible sgoldstino production has been estimated using 1 .5 M generated events; 
in this sample the sgoldstino mass mp is varied from to 200 MeV/c^ in steps of 10 MeV/c^. These 
simulated signal events were weighted using the matrix element given in [2|. 

Some information on the data processing and 
reconstruction procedures is given in E [HI |9||, 
here we briefly mention the details relevant for 
7r~7r° -I- missing energy event selection. ,^6 

The muon identification (see [8|) is based on 
the information from the SPi and the HC. The elec- 
tron identification (see [7|) is done using the E/p o 
ratio — the energy of the shower associated with a 2 /o' 
charged track and the charged track momentum. | 

A set of cuts is developed to suppress back- u 
grounds to possible sgoldstino production: 

0) Events with one reconstructed charged track 
and with two showers in the calorimeter SPi are 
selected. We require the effective mass m(77) to w 
be within ±50 MeV/c^ from m^o. Events with a 
vertex inside 400 < z < 1650 cm are selected. ; 

1) "Soft" charged pion identification is applied: 
tracks identified as electron or muon (as described 
in QIHl) are rejected. 

2) Events with missing energy = Eteam- Figure 3: Missing energy spectra, second run 
E„- - E^o above 3 GeV are selected. This cut data. Points show real data, upper histogram — 
serves to reduce the K„2 contamination. The miss- background MC events, lower histogram — sig- 
ing energy spectra for the second run data, MC nal MC events for mp = 90 MeV/c^. 
background and MC signal with m p =90 MeV/c^ 

are compared in Fig. 3. 

3) To reduce the background caused by secondary charged particles in the photon sample, the MH 
information is used. The distribution of the distance Ar between an SPi shower and the nearest MH 
hit in the plane transverse to the beam is shown in Fig. 4. An event is selected if at least one shower 
has Ar greater than 10 cm. 

4) The events where one of the photons is suspected to be irradiated by the charged particle in 
a detector material upstream/inside M2-magnet are rejected. Such photons have nearly the same x 
coordinate (in the direction of the magnetic field) as the charged track on the SPi, i.e. an event is 
rejected if for one or both showers have the difference Ax = \xch — x-y\ < 7 cm . 

5) Events which pass Kes 2C-fit Hi are removed. 

6) An additional cut on the charged pion momentum p* in the kaon rest frame (vr~ ) < 180 MeV/c 
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is applied to suppress the tails from the ir^2 decay. 

7) The pion identification, mainly based on the HC 
information is required . The efficiency and muon sup- 
pression of this selection has been determined using 
Kij,2 and Kt^2 decays. The efficiency is found to be 
70% and 80% for the first and second run data respec- 
tively. The remaining muon admixture is equal to 3% 
and 2%. 

8) The events with missing momentum pointing to 
the SPi working aperture are selected in order to sup- 
press 7r~7r°7 background. 

9) The "Veto" cut uses information from the Guard 
System (GS) and the guard electromagnetic calorime- 
ter (SP2): absence of signals above noise threshold is 
required. 

The data reduction information for the second run 
is given in Table 1 and is compared with the MC- 
background and with the MC- signal for the sgoldstino 
mass of mp = 90 MeV/c^. The cut suppression fac- 
tors Wi-i/wi for the MC data are calculated using the 
corresponding matrix elements. 

The influence of the last cuts on the missing mass 
squared spectra (Pk — P^r- — P-k°Y is shown in Fig. 5. The left wide bump in Fig. 5 is due to 
Kn2, decays, the shift to the negative missing mass squared is caused by the use of the pion mass in 
its calculations. The second peak is caused by 7r~7r°(7r°) decays with photons from the second 7r° 
escaping detection both in the SPi and in the "Veto" system. 




40 60 
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Figure 4: The distance Ar between the SPi 
shower and the nearest MH hit, second run 
data. 



3 Analysis and results 

As a result of the previous cuts, especially the last "Veto" cut, the 7r°-signal practically disappears 
from the 77 mass spectrum for certain missing mass (mp) intervals. To illustrate that, the m(77) 
spectrum after cut (2) is shown in Fig. 6 (left half). The right half of Fig. 6 shows the 77 mass spectra 
after the "Veto" cut for several mp intervals. It is clearly seen that the 7r° signal survives for the region 
of the 7r~7r°(7r°) decay only. This situation allows us to perform an effective background subtraction 
procedure: the 771(77) spectra for 10 MeV/c^ mp intervals, starting from mp = 0, are fitted by the 
sum of a Gaussian (with the fixed width of 9.3 MeV/c^ and the fixed mass of 135 MeV/c^) and a 
quadratic (or linear) polynomial. Such fits are shown in Fig. 6 (right half). Fig.7 shows the obtained 
numbers of 7r° (N^o) plotted in the corresponding missing mass (mp) bins. A possible constraint 
A^^o > is not imposed during the fit, thus A^o is negative for some bins. This was done to avoid 
a specialized treatment of asymmetric errors during later analysis. The sgoldstino signal would look 
like a peak in this "filtered" missing mass spectrum of Fig.7, similar to that of a well seen 7r° peak (at 
135 MeV/c^) from the K' 7r-7r°(7r°) events which pass the "Veto" cut. 
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Table 1 : Event reduction statistics for the 2^°- run, the data, the background MC and the signal MC 
with nip = 90 MeV/c^. 



Cut 


data 




BG MC 


"',-1/ 


signal MC 


"',-1/"'/ 


0) 1 Ti^, near ///(tt") 


9943046 




5512890 




98289 




1) no (e, pL) 


7771606 


1.28 


4545059 


1.19 


93470 


1.05 


2) E^i, > 3.0 GeV 


1123220 


6.92 


588735 


7.78 


82602 


1.17 


3) MH filter 


939052 


1.20 


516922 


1.19 


74744 


1.09 


4) conv. gammas 


722622 


1.30 


426286 


1.22 


56513 


1.25 


5) no fit 


458338 


1.58 


201580 


2.27 


35906 


1.69 


6)p*{n-) < 180MeV/c 


326935 


1.40 


134706 


1.28 


35698 


1.01 


7) 7r~ identification 


122804 


2.66 


68380 


2.06 


33401 


1.06 


8) 10 < r < 60 cm 


108992 


1.13 


60698 


1.12 


31431 


1.06 


9) Veto 


31451 


3.47 


18674 


3.58 


31104 


1.01 



data MC BG 




Figure 5: Missing mass squared distributions after cuts 6 (a) , 7 (b) and 9 (c), left column — sec- 
ond run data, right column — MC for dominant K~ decay modes. The bin size is equal to 0.0005 
(GeV/c2)2. 
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Figure 6: a) m^'yy) spectrum for events with Emis > 3.0 GeV with the result of the fit by the sum of 
the Gaussian and a third degree polynomial; ^(77) spectra after the "Veto" cut for the missing mass 
mp intervals 120-130(b), 130-140(c), 140-150(d) and 150-160 MeV/c2(e). 



No significant signal is observed. To obtain the upper limits for the number of sgoldstino events at 
different mp (Ngig), the "filtered" missing mass spectrum is fitted by the sum of the signal Gaussian 
(with a fixed width of a = 11.1 MeV/c^, as determined from the signal MC) and a background. The 
background is described by the sum of two components: a Gaussian for the 7r^7r°(7r°) peak plus the 
constant term. The integral of the signal+background function over a particular bin (-i-) enters into 
the function which is minimized by the "MINUIT" program [11|. Fig. 7 shows an example of 
such a fit with mp = 185 MeV/c^. To be consistent, we do not impose the constraint Nsiy > 0. We 
checked that the constrained fit gives comparable or smaller error estimates (corresponding to lower 
upper limits), i.e. the procedure used is more "conservative". One-sided upper limits for the number 
of signal events at the 90% confidence level are calculated as 

NuL = max{N,,g,0) + 1.28 -a , (3) 

where a is an error estimate for Ngig. The upper limit for the branching ratio (UL) is calculated as 

^ NuL- 0.2116 ■eiK^2) 

where 0.2116 is the branching ratio Br{KT,2), and N{Kt,2) is the number of reconstructed 
7r~7r° decays found to be ~ 1.5 M events for the first and ~ 4.5 M events for the second run. The val- 
ues e and e{Kn2) are respective efficiencies for n^n° P and Kt,2 decays, which include both 
the reconstruction efficiency and geometrical acceptance. To avoid the systematics associated with the 
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"Veto" cut (cut 9), we have applied it to the selected Kt,2 events. The measured efficiency of this cut 
is 90 % and is explained by random signals in the Veto system due to the presence of the beam halo and 
event overlaps. 

Since the Veto inefficiency has the same influ- 
ence on both the signal and the normalization 
sample, it cancels out in (4). The signal ef- 
ficiency e rises monotonically from ~ 3.3% 
to ~ 11.3% in the mass interval from to 
170 MeV/c^ and then drops slightly to the 
value ~ 10.3% at 200 MeV/c^. This behav- 
ior is practically identical for both runs. The 
Kt^2 efficiency e{KT,2) is equal to 15.8 % for 
the first run and 22.4 % for the second run. 
The weighted average ratio ejeiK-^^) with the 
weights proportional to the run statistics is 
used for the combined data sample. 

The final results obtained with the com- 
bined statistics of two runs are given in Fig. 8. 
The left part of this figure shows a compari- 
son of our results with that published by the 
E787 coUaboration Q. 

The obtained upper limits can be trans- 
formed into the limits for the value of the 
modulus 



is 300 
(U 

> 



200 — 



100 — 



I '''12 I 



(see equations 1,2). The corre- 
sponding limits are compared in Fig. 8 (right 
part) with the theoretical limit evaluated using 
Kl — Kg mass difference. A generalization 
of the formula B9 in the Appendices B of the 
paper [2| for the non-zero nip |12| was used 
when extracting the | /if^ | . 

The systematics in the upper limits is dis- 
cussed in the following section. 




missing mass m(P), GeV/c 



Figure 7: The number of n° from the fit of the 7/1(77) 
mass spectra plotted in the corresponding missing 
mass nip bins. The curve is the result of the fit with 
an assumed sgoldstino signal at nip = 185 MeV/c^ . 



4 Systematics 

The comparison of our data with the background MC shows some discrepancies: although the total 
suppression factors of the selection cuts for the data and MC are equal within a few percent, a signifi- 
cant difference of 25-30 % is observed for cuts 6,7(see Table. 1). Some differences seen in Figs 3 and 5 
are explained by the non-Gaussian tails (absent in the MC data) of the beam momentum distribution 
and higher noise levels in the real data. 

To estimate the systematics related to the observed differences, we have calculated from our data 
the well known branching ratio Br{K~ n^n°n°), using the events of this decay with one n° 
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Figure 8: The 90 % CL upper limit for the Br{K^ ■n^-n°P) versus sgoldstino mass compared 
with the E787 upper limit(left), the 90 % CL upper limit for the \h^2\ compared with the theoretical 
limit from — Kg mass difference(right). 



missed, i.e. events with the topology K~ 7r~7r°(7r°). In this test the missing 7r° plays the role of 
sgoldstino. 

For the test, we have repeated the every step of our analysis, i.e. event selection, which includes: 
cuts 0-8, filtering of the missing mass spectrum, fitting and normalizing to the — > ii^tx° decay. 
The only cut which was removed is the "Veto" cut (cut 9) in order not to suppress the escaping n° . 
The branching ratio Br{K~ 7r~7r°7r°) has been measured to be 1.64 ± 0.03%(stat.) compared with 
the PDG value 1.73 ± 0.04 % lfT51 . i.e. the observed systematics is ~ 5%. It has been taken into 
account by the multiplication of the UL in (4) by the factor of 1 .05 for all data points. 

This procedure ensures the absence of the systematics in the upper limit for sgoldstino mass 
mp ~ m^o . To estimate the residual systematics for other mass bins the upper limits were recalculated 
for: 

a) different cut levels in the cuts 2 — 9; 

b) with a different background parametrization for ^(77) (Fig. 6) and N^^o (Fig. 7) spectra; 

c) variation of the theoretical parameters | 2| used in the signal Monte-Carlo program. 

For each new set of cuts and for each mass point nip the upper limits were calculated exactly in 
the same way as for the main set of cuts. The observed variation of the results gives the estimate of 
the remaining systematics. For example, for mp = 95 MeV/c^ we got 17%, 11% and 4% variation 
for the contributions a), b) and c) respectively. 

This systematics for each point was taken into account by the corresponding increase (by the 
summary multiplicative factor) of the upper limits ( 21% increase for mp = 95 MeV/c^ equal to sum 
of the contributions a), b) and c)). For other data points the systematics varies from 16 to 25 %. 
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5 Summary and conclusions 



A search for the light pseudoscalar sgoldstino in the -k^-k°P decay is performed. It is assumed 

that sgoldstinos decay outside the detector. No signal is seen in the mp mass interval between and 
200 MeV/c^. The upper limits at 90% CL are set to be ~ 9.0 ■ 10^^ for the sgoldstino mass range 
from to 200 MeV/c^, excluding the interval near m{Tc°), where the limit is ~ 3.5 • 10^^. These 
results improve the confidence limits published by the E787 collaboration. Our limits improve also 
the theoretical constraints on from Kl — Ks mass difference. 

The authors would like to thank D.S. Gorbunov, V.A. Matveev and V.A. Rubakov, for numerous 
discussions. 
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